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Abstract. Perturbation theory within Newtonian approximation is presented for 
cosmological models with varying physical constants. Analytical solutions for 
perturbations dynamics are obtained for each Gurzadyan-Xue model with pressureless 
matter and radiation. We found that perturbations grow during entire expansion 
within GX models, including curvature- and vacuum-dominated stage when they cease 
to grow in the standard cosmological model. 



1. Introduction 



The formula for dark energy, derived by Gurzadyan and Xue predicts the observed 
value for the density parameter of the dark energy without any free parameters [I] . The 
formula reads: 

71 ?IC 1 7T C 4 1 

PGX = 8Lf^ = 8G^' ( ' 

where h is the Planck constant, the Planck length is L p = (hG) 1 c~ 3 / 2 and c is the speed 
of light, G is the gravitational constant. Here a is the upper cutoff in computation of 
vacuum fluctuations and we take it to be the scale factor of the Universe, although other 
choices are possible L2J. According to Zeldovich the vacuum energy (0) corresponds 
to the cosmological term 

= **(;)'. 0) 



In General Relativity all the quantities in ((TJ and (J2J) are constants, except for the 
scale factor, which is a function of cosmic time. Therefore, adopting Gurzadyan-Xue 
(GX) scaling (JTJ) one has to consider possible variation of physical constants||. 

| e-mail address: veresh@icra.it 
§ e-mail address: gegham@icra.it 

|| Generally speaking, the vacuum energy is known to be constant only in Minkowski spacetime. 
However, assuming it is constant one has unique connection between the vacuum energy density and 
the cosmological term and can interpret one of them as a fundamental constant of physics. 
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Simple cosmological models, based on above ideas we proposed in |4j. In particular, 
models with varying cosmological term, speed of light and gravitational constant were 
considered. Qualitative analysis of solutions in these models revealed some interesting 
features, in particular the presence of separatrix in the phase space of solutions. This 
separatrix divides the space of solutions into two classes: Friedmannian-like with initial 
singularity and non-Friedmannian solutions which begin with nonzero scale factor and 
vanishing matter density. Each solution is characterized by a single quantity, the density 
parameter which is defined in the same way as in the standard cosmological model: 

J L m i 
He 

where fx is matter density, H is Hubble parameter, fx c = is critical matter 
density, and index "0" refers to the values today. Separatrix in all models is given by 
the density parameter Q m ^2/3 that indicated some hidden symmetries between the 
models, although cosmological equations look very different. Later in JT] the origin of 
these symmetries was discovered and attributed to existence of invariants in GX models. 

Note, that the energy-momentum conservation does not lead to mass conservation 
of the usual matter in models with varying constants. This can be interpreted as creation 
of matter due to variation of constants. As a consequence, the usual power law fx oc a -3 
does not hold in expanding Universe. In general fx(a) is obtained by solving cosmological 
equations. 

Analytical solutions for all GX models both for matter density and the scale factor 
are obtained in j^j. It turns out that the most simple solutions for the scale factor 
are again those of separatrix. In one model it is exponential, in the others they are 
polynomials. 

Simple GX models were generalized to include radiation in [G], where analytical 
solutions were also obtained. The main difference in that case is absence of separatrix 
for matter density solutions in models with radiation: all solutions (with exception 
of model with varying cosmological constant) begin with vanishing matter density in 
contrast to Friedmannian solutions. There is, however, a particular density parameter 
in those models that divides solution again into two classes: in the one class the matter 
density increases in matter-dominated epoch, in the other it decreases. 

Thus, GX models possess interesting and nontrivial features. The original 
motivation for these cosmological models, however, is the fact that the predicted dark 
energy density is close the value inferred from the set of current observations. Along 
this line, it is of crucial importance to further confront predictions of these models with 
observations. We have performed likelihood analysis of supernovae and radio galaxies 
data within GX models in |7 • Such important characteristics of models such as age and 
deceleration parameter were also computed. Our results show that models with density 
parameter smaller than the separatrix value do not pass observational constraints. This 
indicates the preference of high matter densities: GX models with 2/3 < Q m < 1 
pass these simple but important constraints. All these results show that background 
dynamics of GX models is viable in front of recent observations. 
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In this paper we turn to analysis of perturbations dynamics within GX models. It 
is well known that the most stringent constraints on cosmological parameters within 
ACDM models were obtained from observations of large scale structure distribution 
of galaxies and clusters, as well as from cosmic microwave background radiation 
anisotropies. We will consider simple Newtonian perturbations, i.e. those with 
wavelengths well inside the horizon. This allows, nevertheless, to understand the main 
differences with respect to dynamics of perturbations in the Friedmannian models. 

The paper is organized as follows. In sec. 2 GX cosmological models are reviewed. 
In sec. 3 we present cosmological perturbation theory with varying constants in 
Newtonian approximation. In sec. 4 this theory in applied to GX models. In sec. 
5 discussion and conclusions follow. 

2. Gurzadyan-Xue cosmological models 

In this section we remind briefly the content of GX models, for detailed discussion see 



Adopting the scaling (0) and postulating the validity of Einstein equations when 
speed of light and gravitational constant vary with time we obtain four simplest cases 
when the speed of light, the gravitational constant or vacuum energy density vary with 
time in such a way that is satisfied. Below we provide cosmological equations for 
each model. The content of the Universe is supposed to be a pressureless fluid with 
mass density /i and energy density fie 2 as well as a radiation field with energy density 
p r and pressure p r /3. 

We also assume that evolution of radiation energy density is described by the 
following equation 



which means there is no energy flows into radiation field. This assumption allows 
to avoid problems with anisotropies of cosmic microwave background and large scale 
structure data for models with varying constants [ELIJ- 
AH GX model with radiation are characterized by two density parameters: fl m and 



2.1. Model 1, with radiation 

We suppose that neither the speed of light, nor the gravitational constant change with 
time, but the vacuum energy density does, so pax oc a ~ 2 - This is the unique case when 
the radiation does not couple to matter jB]. 

This model is defined by the following cosmological equations |Hj 



[6j. 






(3) 
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Figure 1. Total energy density as a function of scale factor for model 1 with radiation. 
Here ao is scale factor today. Critical energy density is defined as p c — fi c c 2 . In this 
and subsequent figures black curve correspond to density parameter Sl m = 2/3, dotted 
curve is for Q m = 0.3 while firm curve is for il m = 1. Also for all figures k = 0. 

f + 3H " = w(z)' a ' (4 > 

where p r o is radiation energy density today, and effective curvature is k! = k — 
k = —1,0,1. The first equation looks identical to Friedmann equation without 
cosmological constant. The second equation corresponds to continuity equation with 
the source term due to vacuum energy variation. This term is positive during expansion 
and may be interpreted as creation of matter. 

Solutions of these equations are represented at fig. 12 . 1| 13 Matter density and 
energy density of pressureless fluid are related through constant speed of light. 

This model contains radiation-dominated (RD) stage which begins from 
Friedmannian singularity if Q m < Q sep , where density parameter for separatrix is 
^sep = 1 x-k/n* ' or W ^ n finite value of scale factor and vanishing energy density if 
Q m > Q sep - Then RD stage is followed by matter-dominated (MD) stage. On last (kD) 
stage curvature dominates^. 

2.2. Model 2, with radiation 

In this model the speed of light changes with time in such a way that vacuum energy 
density (0) together with (0) remain constants, so that the latter may be interpreted as 
cosmological term. The speed of light increases during cosmic expansion: c oc a 2 . Our 
cosmological equations are 

More precisely, terms of the order a~ 2 dominate in the first cosmological equation. 
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Figure 2. Scale factor as a function of time for model 1 with radiation. Here time is 
measured in inverse Hubble parameter today, for which we take Ho = 70 km/(s Mpc), 
so Hq 1 = 14 GYr for this and all subsequent figures. 
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where A' = A (1 — 3k/ir 2 ) is effective cosmological constant. 
Solutions of these equations are shown at fig. El HI and El 

The difference with model 1 is that now energy density and matter density behave 
differently (see fig. El and HJ since in the simple relation p = pc 2 the speed of light is no 
more a constant but a function of the scale factor, so now p oc pa 2 . Separatrix solution is 
very special since the total energy density for this solution begins in infinity, decreases, 
then vanishes at some value of the scale factor, and increases without bound again for 
large a. For smaller density parameter, than the separatrix one, solutions admit negative 
energy density for finite interval of scale factors, and for larger density parameter the 
energy density is always positive definite. In contrast (as can be seen at fig. EJ), matter 
density is always negative for small scale factors (RD stage). For separatrix solution 
the matter density stays constant during "vacuum-dominated" (VD) (dominance of the 
vacuum energy density in cosmological equations) epoch, which follows after MD stage; 
there matter density decreases for Q m > Q sep and increases otherwise. 

2.3. Model 3, with radiation 

Here gravitational constant decreases with cosmic expansion G oc a~ 2 and it is assumed 
that vacuum matter density pax is a new constant. Cosmological equations read in this 
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Figure 3. Total energy density as a function of scale factor for model 2 with radiation. 
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Figure 4. Matter density of pressureless component in units of critical density as a 
function of scale factor for model 2 with radiation. For solution with Cl m > il sep the 
matter density has a maximum which is not shown here. 
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At fig. 6, 7 and 8 we illustrate solutions of these equations. 
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Figure 5. Scale factor for model 2 with radiation as a function of time. 
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Figure 6. Total energy density as a function of scale factor for model 3 with radiation. 

One can see that during the VD stage both total energy density and matter 
density of pressureless fluid stay constants for separatrix solution; they both increase 
for Q m > Q sep and decrease otherwise. 

2.4- Model 4, with radiation 

In this model it is assumed that the speed of light changes with time, but the vacuum 
energy density pax stays constant. This leads to c oc a 1 / 2 dependence. 
Cosmological equations are 
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where /?=^(Gp GX ) 1/2 (l-f§). 

Solutions in this case are represented at fig. 9, 10 and 11. 

The role of separatrix in this model is similar to the model 3 with respect to the 
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Figure 9. Total energy density in the model 4 with radiation. 
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Figure 10. Matter density in the model 4 with radiation. 

energy density, but the matter density of pressureless matter behaves differently since 
now p oc a/i. 



3. Cosmological perturbations with varying constants in Newtonian 
approximation 

Since cosmological equations and their solutions for GX models are known, it is of 
interest to consider evolution of perturbations in these models. This is relevant for 
understanding the dynamics of large scale structure, its analogy and difference with the 
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Figure 11. Scale factor in the model 4 with radiation. 



standard cosmological model. Due to some peculiarities of GX models this dynamics 
turns out to be quite different. 

We remind that Newtonian equations for perturbations may be used to describe 
evolution of perturbations in nonrelativistic component and for modes well inside 
horizon. 

In general, when both speed of light and gravitational constant change with time 
we have (see Appendix): 
continuity equation 

Euler equation 

dvq dvq dp chp_ _ c . . 

dt Wp 8x? ndx a dx<* " Va c' 1 } 

and Poisson equation 

Ay? - AirGfi = 0, (13) 

where \i and p are density and pressure of the perfect fluid, which is moving with 
velocity v Q in a gravitational potential ip. Both the continuity and Euler equations 
contain additional terms: A* — 2^ in continuity equation and — term in Euler 
equation. This means, that mass in comoving frame is not conserved in models with 
varying constants, so there is continuous creation of matter. Also velocity flow changes 
not only due to pressure and gravitational field gradients and universal expansion, but 
also due to variation of the speed of light. 

Density perturbations in expanding Universe obey second order differential 
equation. In Newtonian approximation, neglecting the sound velocity this equation 
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reads (see Appendix) 

c"/ x c'/ x / 3 1 c if'\ ArrGft ^ 

where 5 = (/i — ft) /ft is dimensionless density contrast, is unperturbed density. 
Here again, in contrast with the corresponding equation for Friedmannian models, 
additional term appears due to variation of the speed of light. These equations are 
well defined to describe perturbations in pressureless component which represents e.g. 
cold dark matter, or baryons after recombination in large wavelength limit. 



4. Linear Newtonian perturbations in GX models 

4-1. Friedmannian models 

Solutions for density perturbations in Friedmann models are well known. During RD 
stage density contrast of dark matter grows + with the scale factor logarithmically*, 

5 oc In a. 

During MD epoch perturbations both in baryonic and dark matter components 
grow linearly with the scale factor, 

5 oc a. 



Finally, in VD, or curvature-dominated (kD) epochs perturbations cease to grow. 
Now we turn to perturbations dynamics specifically within GX models. 



4-2. Model 1 with radiation 

The solution for matter density is 



c 2 tt / a \ f a\ 3 

and one recovers Friedmannian behavior ft oc a~ 3 in the limit of small scale factors. 
However, it is during RD stage where dynamics is determined by radiation. According 
to this solution at late expansion dominates the curvature. 

In this model c = const, G = const, so the perturbation equation is 

-^=»- (14) 

+ In this paper we consider only perturbations with wavelength larger than the corresponding Jeans 
wavelength. 

* We are interested in growing modes only. 
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4-2.1. Radiation- dominated epoch. Within Newtonian approximation one can still 
follow dynamics of perturbations in the pressureless component during RD stage |lUj . 
This results in the following solution 

5 oc In a, 

so perturbations evolve exactly as they do in Friedmannian models. It is not surprising, 
since during this stage we have usual dependence of the matter density on the scale 
factor, and ()14|) is also the same as in Friedmannian models. 

4-2.2. Matter-dominated epoch. Here we take into account only a~ 3 terms in (|4.2jl and 
also put p r o — > 0, k = when calculate H(a). The solution is 

5 oc a, (15) 

and it again coincides with the corresponding solution in Friedmannian models. The 
reason is the same as for RD stage. 

4-2.3. Curvature dominated epoch. Here we account only for a~ 2 terms in (|4.2jl and 
also put p r o — > when calculate the Hubble rate. 
Thus we get the solution 



5ocaHV 5+ ^ V (16) 

Interesting, that for zero spatial curvature this is "golden section" solution 

5 oc aH^-i) (17) 

Perturbations grow at kD stage due to existence of energy transfer from dark energy 
into pressureless matter which manifests itself in the presence of the source in continuity 
equation of this model. 

As we see, in contrast with the Friedmann models, perturbations grow during entire 
cosmological expansion. 

4-3. Model 2, with radiation 

Matter density depends on the scale factor as follows 

/ ti n \ " 

(18) 



A / a \ a 4vr 2 pro 

A* = 1-77 1 — H — ^0 + 



AttG \ a J a 5aoA 



a J 



As in this model c oc a, the perturbation equation becomes 

+ (1 + 1) -^' = °- (19) 

In all GX models with radiation, except for model 1, matter density is negative in 
the beginning of cosmological equations (see figures above). This leads to anomalous 
behavior of perturbations during RD stage and in what follows we do not consider 
dynamics of perturbations for RD epoch. 
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4-3.1. Matter-dominated epoch. In contrast to model 1, matter-dominated epoch 
corresponds to the prevalence of a -1 terms in (fTHj) . We also put A = 0, p r0 — > 
in the background equations when calculate the Hubble rate. 
The solution is 

S oc y/a. (20) 

Thus perturbations growth is somewhat slower comparing to Friedmannian case in 
MD epoch. This is a consequence of matter creation in all space due to variation of the 
speed of light with cosmic time. 

4-3.2. Vacuum- dominated epoch. Here we take into account only a terms in (|18|) and 
also put p r o — > when calculate the Hubble rate. 
In this way we obtain the solution 

5ocaHV^^ V. (21) 

Matter is also created during VD stage but this leads to an opposite effect, 
comparing to the case discussed above, namely it stimulates the growth of perturbations, 
like at curvature-dominated stage in the model 1. 

For zero spatial curvature k — it follows 

5(a) = cia^-i. (22) 
Again, perturbations grow at all stages of expansion. 

4-4- Model 3, with radiation 

Cosmological equations lead to the following solution for matter density 

H = 2p GX + * Oio - 2pox) + ^ f 1 " (^)l • (23) 
a 6c z a \ V a / J 

In this model c = const, G oc a~ 2 , so the perturbation equation is the same as in 
model 1. 

Solutions for perturbations within this model during MD and VD epochs are exactly 
the same as in model 1 during MD and kD epochs respectively. 

4-5. Model 4, with radiation 

Matter density depends on the scale factor as follows 



In this model c oc a^, G = const, so the perturbation equation is modified 
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4-5.1. Matter- dominated epoch. Considering only a 2 terms in (|24|) and assuming in 
addition /3 — 0, p r0 — > in the cosmological equations we find the solution 

b (X a * 4. (26) 

This is slower growth than in the Friedmannian model and model 1, but faster than 
in the model 2. 



4-5.2. (3 -dominated epoch. Taking into account only a 1 terms in (}2"4*|) and also putting 
Pro — > to calculate the Hubble rate we find the solution 



5 oc a V . (27) 

For zero spatial curvature we obtain 

6(a) = a^~\ (28) 
Again, perturbations grow during entire expansion. 



5. Conclusions 



Dynamics of perturbations in Gurzadyan-Xue cosmological models is studied within 
Newtonian approximation. With this goal relativistic hydrodynamic equations for 
models with varying constants are derived. Non-relativistic limit for these equations 
is also obtained, which provided the basis for analysis of perturbations dynamics 
within Newtonian approximation. Analytic solutions are presented for the model with 
pressureless ideal fluid and radiation. In contrast with the standard cosmological 
model perturbations grow during entire cosmological expansion, including curvature- 
dominated and vacuum-dominated epochs. 

Only Poisson equation preserves its form comparing to the Newtonian 
hydrodynamics. Both continuity and Euler equations contain additional terms due 
to variation of the physical constants, namely the speed of light and the gravitational 
constant. Thus both at cosmological and at local level continuity equation contains 
source terms which may be interpreted as creation of matter. 

We found no difference in dynamics of perturbations between our model 1 (varying 
vacuum energy density) and standard cosmological model for radiation- and matter- 
dominated stages. However, on the last stage perturbations grow according to the 
"golden-section" rule. In other models we have found polynomial dependence of the 
density contrast on the scale factor. 
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7. Appendix 

7.1. Hydrodynamic equations 

It is our goal in this section to derive hydrodynamic equations for models with varying 
speed of light and gravitational constant in Newtonian limit. We start with the interval 



„2 



ds 2 = (c 2 + 2ip) dt 2 - dr\ 
where r a = (x, y, z) are cartesian coordinates, so metric coefficients are 

fl'oo = c 2 + 2<p, g u = g 22 = #33 = -1- 
In Newtonian limit 

c 2 

is a small parameter, and in Taylor expansion we will keep only zeroth order terms 
with respect to it. 

Christoffel symbol has the following nonvanishing components 



1 00 



of) 



CC + (p 

c 2 + 2(^' 



■pa 
1 00 



dip 
dx a ' 



dip 



c 2 + 2ip dx a 

Here Latin indexes take values from to 3, while Greek indexes run from 1 to 3; a 
dot denotes differentiation with respect to time. 

Ricci tensor has the following nonvanishing components 



R 



1 



Rn 



1 



d 2 ip 1 / dip 

d(x a ) 2 ~ c 2 + 2ip [dx* 







c 2 + 2<p 



Aip 



1 

c 2 + 2ip 



{Vipf 



;Aip. 
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7.1.1. Poisson equation. From Einstein's equations 



where 

8ttG 



K = 



C 4 



for the fluid at rest one has 
4ttG 

so we obtain the usual Poisson equation 



Rq - — —fi, 



A(p = AirGfx. 
Energy-momentum tensor for a perfect fluid is 
Ti k = cou k Ui-pgi k , 

where u = fic 2 +p, u { = -^— (l,v a ), u { = -^= (c 2 + 2ip,-v a ), 7 = (l - J) 
In components 

T = u-f 2 - p, 

T a = cu 7 2 v a = ~c 2 T a °, 

T Q /3 = -cj(-) v^V/3 -p5 aP . 



. c, 

Euler and continuity equations both follow from energy- momentum conservation 



k ^ \/— g /t dx k \ ' 2 cte* 

We again expand all terms in this equation in Taylor series with respect to small 



parameters 



c 
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7.1.2. Continuity equation. From zeroth component of energy- momentum conservation 
we find 



1 d^/=g~ k 
--i-o 



'—g dx k 

1 {dy/=g 



u rp 



dx k 



2 dt 



r o T o 
_i oo-'o 



_i 00 J o 

d 



dt 

^0 

00 A o 



/V <)l ■ ' <).!■' 

ilp w 7 2 v Q + ( ^ - 4^ J (^ 7 2 - p) + ^ (cu 7 2 - p) + 



<9a; a 

a0^0 



r0 rp , "pO a | ^rp 

1 nnin + 1 „n J n J 



c 2 + 2(p dx c 



+9^ Kv«) 



1 <%> 2 

^7 v a + 



G 

'G 



c 2 <9x a 



9t 
di 



(W p) • Try ' - l> 



■) + 



Taking Newtonian limit of this expression, and dividing all terms by c 2 we find 

'G 
G 



1 dtp 1 (G A c\ 2 1 9 ( o\ 1 d ( 2 \ 



c 2 <9a; a 
c 2 cte' 



'G c\ c da d . . 

--4-U + 2 M - + - + s? (^) 



<9u <9 , . /G „c\ 1 9ip 

SF + ^ ( " V ° )+ (G- 2 cJ" + ^'' V ° = a 

The last term is small comparing to the second, as can be seen from the following 



2y? < c 2 , 
2/iv a (p < /iv a c 2 , 



9 



,9(//V Q ) 



2/iv a — — + 2^ 



< c 



; <9(/xv c 



&r a <9x a <9a; a 

Therefore, continuity equation for models with varying constants in Newtonian 
approximation reads 



d/i d 
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7.1.3. Euler equation. For the remaining 3 components of the energy- momentum 
tensor we have 

1 d^f=g~ o _|_ / rp k\_ ld 9ki T u _ 



-g dx k 



1 ( d^g- dV^g p\ k o + 9_( T o\ , 
a + dx? a J k + dt" a ) + 



d 



-9 



dt 



( rp 0\ _ 1 ^5*00 rpQQ _ pO rp OipO T P i *2_rp Q , ^ ( ' rp 

dxP \ ) 2dx" ~ 00 a + ^ a « a + dt" a 



+ 



1 dip 



CC + (f 
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c 2 + 2<p \c 
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c 2 + 2ip dx@ 
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dt u 
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G 



c 2 + 2ip dx° 
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(^7 2 - p) 
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7 



" ( c ) VaV/? + P6a(s + ^dx^ ^ - 2p ) =0 



Taking Newtonian limit of this expression we get 



c dip VqV/j (G c\ d d 

-Z^a ~ ~ \g~ cl fIVa ~ di (/iVa) ~ dx~P (/iVaV/3) ~ 

1-1)- 



dp dip 



d 



dx a dx a 
dip 



(fJLV a Vp) - 



dp d (/xv Q ) 



dx a 

dv a 
/xv r 



— /XV 



dp dip 

— V^r— = 0. 



c dv Q 
+ /xv Q - - /X— - 
c dt 



dx? dx a ^dx a 

Finally, dividing this by /x we find the Euler equation for models with varying 
constants in Newtonian approximation 

dvq dv a dp dip _ c 
dt V/3 ^ fidx a dx a ~ Wa c 

7.2. Perturbation equations 

First, introduce co moving coordinates, 

r a ax a , 
then for velocity we have 

v Q = Hr a + u a . 
Perturbation equations now read 

% + SHfi + -d a (fiu a ) +fl (<l-2-)= 0, 
dt a \G c I 



(29) 



Dynamics of perturbations in Gurzadyan-Xue cosmological models 19 

cP (x du 111 c 

-rr^a + -T^r + Hu a + -u^dpUa-l d a p+ -d a ip+ (Hax a + u a ) - = 0, (30) 

at 2 at a an a c 

d 2 <p - ^Ga 2 fi = 0, (31) 

where d a = d 2 = Differential operators transform from physical to comoving 

coordinates as > ^d a , — > — Hxpdp. 

Background solution for /i = /l(t) , p = p(t), <£> = and M a = is given by the 
system 



<i 2 a 1 _ cfc 
,,a "I - a a ip tiax a — , 
at z a cat 

d 2 (p - 4:7iGa 2 fi = 0, 
Consider perturbations about the background solution 



11=^(1 + 8), p = p(l + 6p), (p = (p(l + 6(p), 
5 < 1, 5p < 1, dip < 1, <5w Q < 1. 



5Ur 



From (jlU) and (JH2J) we have 



Then 



d 1 

— |/2 (1 + 5y)} + 3H [/2 (1 + 6fi)] + ~a Q {[/2 (1 + <5w Q } 

^^7 + AiTTT + 3 #A^ + [0U a ) = -fub 

at at .a 
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dfx 
~dt 



3Hfi + fi 







"95 




+ /2 


9t 



85 1 n . . , 
— + ~a Q (du a ; 
at a 



0. 



1 



<9 a (5u c 



0. 



(32) 
(33) 



(34) 



d a [p(l + + 



From (J3UJ) and (J32*j) we obtain 

d 2 a d5u a TTC 1 

+ + Hdu a + 

at 2 at a/i(l + d)) 

1 dc 
+-d a [<p(l + 8(f)] = -— (ifax a + 5u a ) , 

— ^— + Hbu a + —d a p + -d a (dip) = -bu a — . 
at an a cat 

Assume, as usual, the speed of sound is v 2 s = so ^d a (5p) = v 2 d a 8, so 

db~u i) 2 1 dc 

— + H5u a + — d a 5 + -d a {Sip) = -Su a —. 
at a a cat 



(35) 
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From (J31)) and (|32j) we have 

d 2 8<p - 4nGa 2 n5 = 0. 
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(36) 



Combine (pM j) • . (jSBj) to a second order differential equation 



d 2 S d ri 
at 1 at La 

= Trpr H a a ow a + -d, 

ct a a 

d 2 o 



d5u a 



H If v 2 1 dc" 

d a 5v a d a H5u a + —d a 5 + -d a (Sep) + Su a — 

a a \ a a cat / 



dt 2 

^-(2H + -\ ~d a 5u a - ^d 2 5 - AnGfc5. 



Finally, equation for density perturbations is 



dt 2 



+ 2H + 



cj dt 



- -^d 2 5 - AirGfLd = 0. 



The same equation may be rewritten where all functions depend on scale factor 
instead of time with conversion of the derivatives 



5 
6 



Had (a) 
H 2 a 2 5"(a) + 



H (a)Ha 2 + H 2 a 



5(a) 



Then density perturbations satisfy the following equation 



5 (a) +5 (a) 



'3 H 1 c'(a) 
a H Ha c 



H 2 a 



-d 2 S- 



AirGfj, 



H 2 a 2 



5 = 0. 



(37) 
(38) 



(39) 



Comparing to the same equation in Newtonian cosmology jTU] , we find additional 



term - which is due to variation of constants. 

c 



